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ABSTRACT 


Supported vanadium and chromium oxide on TiOa support were considered to understand 
the effect of surface metal oxide loading and type of metal oxide species on the propane 
ODH reaction. It was observed that for the sub-monolayer and monolayer catalysts, the 
propane activity increases with an increase in loading of surface metal oxide sites. The 
propene selectivity, however, followed different trends. To understand these two effects: 
the increase in surface metal oxide species and change of surface metal oxide species, the 
kinetic parameters were estimated for several models using Genetic Algorithms (GA). 
Using GA, a potential technique for parameter estimation in non-linear kinetic models, 
kinetic parameters were successfully determined for the different t5q)es of power law 
models. Based on the kinetic parameters the best-suited model for each catalyst is 
selected using statistical criteria based on F-test for simultaneous reactions. It was 
observed that the 3 % V205/Ti02 and 1 % Cr203/Ti02 are best explained by the PL -1 
model for which CO and CO2 are secondary products. However, for 4 % V205/Ti02 and 
4 % Cr203/Ti02 catalysts, PL -4 and PL -3 models are best suited to explain the 
experimental data. The PL -4 model suggests that CO is produced as primary and 
secondary products and CO2 is formed by the oxidation of CO. The PL -3 model 
considers that CO is secondary product and CO2 is a tertiary product. Based on the 
kinetic parameters of these models the effect of loading and type of the metal oxide 
species is successfully explained. The main considerations being the nature of the pre- 
exponential factors and activation energies. 
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Chapter 1 

INTRODUCTION 


Traditional sources for lighter alkenes are mainly dehydrogenation and cracking 
processes. However, these are relatively expensive processes. With the growing demand 
in petrochemical industries for lighter alkene feed stocks, it is imperative that new 
processes are developed. This is the impetus behind the search for potential alternatives 
to replace the existing traditional sources. Alternative processes like lighter alkane 
oxidative dehydrogenation, though still at the stage of infancy in terms of its commercial 
viability, appears to be more promising for various reasons of interest [1-2]. 

Oxidative dehydrogenation (ODH) of alkanes to alkenes is thermodynamically 
favoured due to formation of water, which is a stable product. An added advantage of 
ODH lies in its irreversible nature and, thereby, theoretically one can achieve complete 
conversion (100%) as there are no equilibrium constraints [2]. On the other hand, the 
endothermic nature of traditional processes, which are feasible at high temperatures 
(-1073 K), have to meet high energy requirements in contrast to the exothermic nature of 
ODH reaction, which is possible even at the temperature as low as 633 K. Furthermore, 
presence of oxygen removes deposited carbon on catalyst during the reaction resulting in 
less deactivation. However, there are certain hurdles to be surmounted before being a 
commercially feasible process. The activation of alkanes is more difficult than that of 
corresponding alkenes. Consequently, the alkene formed is more prone to degrade in the 
presence of oxygen to carbon oxides. Thus the problem lies in maximizing the formation 
of alkenes. Ideally, it requires a catalyst that would significantly accelerate the abstraction 
of hydrogen from the hydrocarbon molecule and hinder both the nucleophilic insertion of 
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oxygen into the molecule and the electrophilic attack of oxygen molecules on the C-C 
bonds [3], 

In an attempt to find high selective catalysts, vanadium oxide (vanadia) supported 
catalyst is found to be a well-established catalyst for partial oxidation of hydrocarbons. It 
is widely used in the petrochemical and environmental industries. Existence of two- 
dimensional active species of vanadia below monolayer coverage is primarily responsible 
for the wide variety of applications. These applications include the ODH of propane, 
which is proposed to be activated by bridge V-O-V bridging bonds [4], The catalytic 
activity of these oxides depends on their surface structure, character of chemical bond 
and the coordination unsaturation of the surface atom. The nature of the vanadium oxide 
species in supported vanadium oxide catalysts are known to depend on the loading and 
the support. Various characterization techniques have been applied to study the supported 
vanadium oxide phases on different metal oxides, which include, Raman spectroscopy, 
solid-state ^ V NMR spectroscopy, UV-vis DRS, and EPR studies. Detailed discussion on 
these spectroscopic studies is available in the literature [4-7,29]. 

Several recent studies are reported on ODH of isobutane and butane using 
chromium oxide (chromia) supported catalysts [8-11]. Similar to the supported vanadium 
oxide catalysts, special catalytic properties of chromium oxide species are also observed 
due to the interaction between the support and the chromium oxide species. Furthermore, 
similar to supported vanadium oxide catalysts, molecularly dispersed chromia species are 
more active for the ODH of propane to propene as compared to bulk Cr 203 [36]. 
Complete understanding of structure-activity relationships requires information firom 
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spectroscopic techniques, such as Raman, XRD, XPS, EPR and TPR along with reaction 
studies [32-33], 

The effect of the type of the surface metal oxide phase (vanadium or chromium 
oxide) and the effect of the surface metal oxide loading are the factors that are required 
for the understanding of the ODH reaction on the metal oxide species. 

Previous reaction studies using supported metal oxide catalysts are performed 
fulfilling various objectives, but few studies are devoted to understand the kinetic 
behavior and mechanism of ODH of Propane. To better understand the effect of loading 
and surface metal oxide species, kinetic modeling can be done. Different models include 
both power law typ and mechanistic models. Mechanistic models, such as Mars-van- 
Kxevelend type redox mechanism [12,13], Rideal type mechanism [14] and mechanism 
involving pseudo-first-order approximation [15,16], have been proposed in order to 
explain kinetic behavior of ODH reaction. Detailed modeling of the ODH reaction, 
without pseudo-first-order approximation, contains highly non-linear functions 
containing local minima. Traditional methods require good initial guess in order to get 
global optima solutions. Genetic Algorithms (GA) provides a potential alternative to 
achieve global minima even when the range to be considered for parameters is extremely 
wide and local minima exists in the parameter space [41]. Performance of GA works on 
the principle of natural evolution. It also depends on parameters like crossover and 
mutation probabilities. 

In the present study, kinetic parameters are estimated for a series of titania 
supported vanadium oxide catalysts and a chromium oxide on titania catalyst using rate 
expressions derived based on the Power Law (PL) type models. Kinetic parameters are 
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determined simultaneously solving non-linear mass balance equations and integrating the 
rate equations over entire mass using fourth order Runga-Kutta algorithm. These 
estimated parameters for different PL models are based on the minimization of multiple 
objective fimctions for simultaneous reactions. Parameters estimated in this way will 
reveal some mechanism involved in each catalyst. 

1.1 Thesis objective 

The present study addresses the catalytic activity of two important catalysts, namely, 
vanadium oxide and chromium oxide on titania support. As mentioned above two factors 
considered are the nature of the surface metal oxide phase and loading of the metal oxide 
phase. For both the catalysts, the surface metal oxide species is considered since these 
species form the active phase for several oxidation reactions. Several loadings of 
vanadium oxide are synthesized on the titanium oxide support. Chromium oxide 
supported on titania catalysts are considered to determine the effect of type of the metal 
oxide phase. The reaction parameters are obtained over the titania supported vanadium 
and chromium oxide catalysts using Genetic Algorithm (GA). Models are discriminated 
using statistical F-test criterion for simultaneous reactions [48]. Based on the results, the 
governing factors controlling the pathways of ODH reaction will be achieved. 

2 Thesis organization 

I s thesis contains 5 chapters. In Chapter 2, a brief review of the literature related to the 
sent study is presented. In Chapter 3, details regarding preparation of the catalyst, 
Derimental set up and the calculations are given. Brief discussion on modeling using 
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GA is also reported. The results obtained from the reaction studies of supported vanadia 
and chromia catalysts on the ODH of propane are presented and discussed in Chapter 4. 
In Chapter 5, conclusions and recommendations for future study are provided. 
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Chapter 2 
LITERATURE EEVIEW 


Conversion of lighter alkanes to their corresponding valuable alkenes is an important 
process due to the increasing interest in the direct use of alkenes as potential feedstocks 
Alkenes find a variety of applications in the petrochemical industries. Studies in the past 
have primarily primarily dealt with using cobalt, nickel and antimony - molybdates based 
catalysts for oxidative dehydrogenation (ODH) of butane to butadiene [17-19], Other 
catalysts, such as, noble metals and metal oxides are successfully tried for several 
dehydrogenation/oxidative dehydrogenation reactions. Metal oxides, such as, Mn02, 
V2O5, C 03 O 4 , CuO, Cr203, M 0 O 3 , Ti 02 and Zr02 as bulk, as binary or in multiple 
composition are commonly employed for the oxidation of hydrocarbons [20-21]. Several 
attempts have been made to raise the yield of alkenes for ODH reaction. However, the 
yield of alkenes is not up to commercial exploitation. In these studies, chromium and 
vanadium oxide catalysts are showing promising activity. Brief summary regarding 
structure and activity of vanadium oxide (vanadia) and chromium oxide (chromia) 
catalysts based literature review is presented along with an account on modeling of ODH. 

2. 1 Vanadium oxide based catalysis 

Vanadia based catalysts have been extensively studied to determine their applicability in 
the ODH of propane reaction. Knowledge on structure-reactivity properties is essential 
for improving the performance of the catalyst and obtaining the fundamental information. 

Deo and Wachs [7] have explained the molecular structures of the two 
dimensional vanadium oxides over layers on different oxide supports including titania 
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under ambient conditions using Raman spectroscopy. They have found that the surface 
vanadium oxide molecular structure depends on the net pH at point of zero charge of 
oxide surface. Furthermore, they observed that a strong influence of the oxide support, 
which was due to the difference in the V-O-support bond strength of the dehydrated 
surface vanadia species. 

Wachs and Deo [23] have proposed several fundamental ideas regarding the 
structure and reactivity of vanadium oxide species on the oxide supports. Ideas regarding 
vanadia species include information regarding monolayer surface coverage, stability of 
monolayer, oxidation state, molecular structure, acidity and reactivity. They have also 
described physico-chemical characteristics of vanadia catalyst compared to other oxide 
catalysts. 

Information regarding the structure of vanadia species, which decide the 
activity of ODH on titania support, is veiy essential. Grzybowska-Swierkkosz [24] has 
studied the distribution of various vanadia species within the monolayer coverage. 
Various types of species include: isolated monomers of VOx of oxo or hydroxo type, 
polymeric species, V2O5 amorphous and crystalline forms. Types of vanadia species and 
their distribution depend on the morphology (rutile or anatase) of the titania support. 

Bell et al [25] have concluded that the structure of vanadia dispersed on 
zirconia is strongly dependent on the surface density and calcinations temperature. It was 
also concluded from the analysis of catalytic and spectroscopic data that polyvanadates 
species, which are highly dispersed on zirconia, are responsible for the ODH of propane 

Pietrzyk and Genser [26] studied the propane ODH reaction using 
approximately five theoretical monolayers of V2O5 on titania (anatase) support xmder 
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transient conditions. Regarding the prospective application of V205/Ti02 catalyst for 
ODH of propane reaction in a circulatory bed reactor, these catalysts, besides requiring 
relatively low temperatures, have the advantage of good stability under the conditions of 
repeated oxidation-reduction steps. 

Gao et al [27] have studied the reactivity properties of the surface vanadium 
species by comparing the structure of the surface vanadium oxide species with their 
corresponding reactivity and selectivity information. In their study it was found that the 
V=0 bond does not influence the reactivity properties of the surface vanadia species 
during hydrocarbon oxidation reactions. The bridging V-O-V bond, however, influences 
some reactions, and the V-O-support bond are found to be the most critical bond 
influencing the reactivity properties. 

Reddy [ 28 ] characterized the V2O5 supported on Ti02 catalyst. Based on 
Raman spectroscopy, it is observed that monolayer coverage occurs at 4%V205/Ti02. 
During the propane ODH reaction, the activity of propane increases up to monolayer 
coverage and then decreases. 

Banares and Cortez [29] have applied an advanced in situ Raman 
spectroscopic online measurement analysis to study the activity and structure of 
vanadium oxide during propane ODH reaction. This approach is considered to be a 
reliable since the structure-activity relationship determined during the reaction. This 
report concludes that surface polymeric vanadium oxide species are more reducible than 
isolated surface polymeric species. 
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2. 2 Chromium oxide based catalysis 

Chromia based catalysis finds varieties of applications. Of these supported chromium 
oxide catalysts have been employed in several chemical reactions due to the simultaneous 
existence of various oxidation states. It is observed that the surface chromia species are 
the centers for catalytic activity. A brief summary of literature review on reactivity of 
chromium oxide catalyst is reported below. 

Hardcastle and Wachs [30] studied the interaction of chromia with TiOj, 
AI2O3, and SiOa under ambient conditions for loadings ranging between 0.5 to 5wt % of 
Cr203. They concluded that surface chromium oxide is present as monomers and dimers 
on alumina, as monomers and possibly dimers on titania, and as monomers and polymers 
on silica. Crystals of Cr203 appear in addition to the surface chromia species above 
monolayer loadings. The differences in distribution of various surface species as well as 
the formation of crystals occur due to the varying surface -support interaction. Based on 
the Raman spectroscopy. 

Vuurman et al. [31] studied the different structures of chromia on alumina 
support. Based on the Raman spectroscopy, they reported that after impregnation and 
drying under vacuum at room temperature, chromium oxide is present as chromate for 
l%Cr/Al203, as dichromate between 5 and 15% Cr/Al203, as trichromate for 20% 
Cr/Al203 and as crystalline Cr03 and trichromate at 30% Cr/Al203. 

Sohn et al [32] observed a strong influence of chromium oxide on zirconia 
support. The surface area of CrOx/Zr02 increased with chromia loading. The authors 
attributed this phenomenon to the strong interaction, which prevents the support and 
chromia from sintering. Upon addition of only small amounts of chromia to zirconia, 
both the activity and acid strength significantly increased. Furthermore, Sohn et al. [33] 
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have reported that addition of chromia transforms the Zr02 support from amorphous to 
tetragonal and monoclinic phase due to strong interaction between chromia and zirconia. 
It was concluded by Wachs et al. [ 34 ] that specific oxide support is a critical parameter in 
deciding the reactivity of surface metal oxide species. 

Scharf et al. [ 35 ] reported that on the surface of supported chromia at least 
three oxidations states (Cr(in), Cr(V) and Cr(VI)) coexist and relative contribution of 
three oxidation states change with loading, calcinations treatment and storage under 
atmospheric conditions. It was also reported that prolonged high temperature calcinations 
results in the decrease of chromia related bands in the Raman and DR. spectra and an 
increased intensity of support related vibrations. They infer solid -state reaction between 
chromia and the support occurs which results into the incorporation of chromium into the 
lattice of the oxide support (Ti02). 

Cherian [ 36 ] studied the structure-reactivity properties of supported chromium 
oxide catalysts for propane ODH reaction with respect to support, surface area, modifiers 
and loading. It was observed that monolayer catalysts are the best for ODH of propane in 
terms of activity and selectivity. Furthermore, the monolayer Cr203/Ti02 was more active 
than the Cr203/Nb205 catalyst, which was more active than the Cr203/Al203. 

2.3 Kinetic modeling 

Kinetic modeling provides us fundamental understanding regarding the roles of surface 
metal oxides on oxide support and surface metal oxide loading. Estimation of pre- 
exponential factors within the monolayer coverage provides the information about the 
variation of intrinsic kinetics of the reaction with loading. Furthermore, estimation of 
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activation energy will enable to predict kinetics route of a particular reaction parameters. 
Recently, genetic algorithm is considered an important tool for the estimation of reaction 
. It is considered to be robust and efficient. A summary of literature survey on ODH and 
parameter estimation by genetic algorithm is given below. 

Parameter estimation using traditional methods, such as, Levenberg- 
Marquardt’s methods dependent strongly on initial guess values. Different initial guess 
values may give rise to different solutions and are explained in the ref [37]. On the other 
hand, objective functions based on non-linear models are common features in most recent 
parameter estimation problems. These nonlinear models using experimental data, 
generally contain more than one minimum for which tedious algorithms are required 
[38], These algorithms are sometimes susceptible to non-global optima. Thus, there is a 
need for development of efficient algorithms to replace the traditional methods. One such 
class of algorithms known as Genetic Algorithms (GA) shows promismg results in terms 
of its accuracy and efficiency. Theoretical foundations of genetic algorithms depend on 
binary string representation of the solutions and the notion of schemata, template 
allowing the exploration of similarities among the chromosomes. This schemata 
representation has led to development of the schemata theorem [39]. 

In an attempt to find initial estimates of rate constants for non-linear chemical 
kinetics. Wolf et al. [40] have exploited GA techniques without a priori assumptions of 
rate determining step in order to apply for the wide range of partial pressures of reactants 
and temperatures. Elliot et ut/ [41] have used the inversion procedure of genetic algorithm 
to estimate rate parameters corresponding to product species measurement data from 
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combustion of fuel experiments. This study suggests that its wide application to other 
chemical kinetics and optimization of other higher order kinetics is possible. 

Moros et al. [42] have estimated the kinetic parameters for the 
dehydrodimerization of methane by using genetic algorithm. It is observed that the 
application of genetic algorithm to generate a good initial point for the non-linear local 
convergence method requires less computing time and increase the reliability of the 
parameters 

Polifke et al. [43] have employed genetic algorithm to investigate the kinetic 
parameters of simplified reaction mechanism for methane combustion. It is also 
determined that using genetic algorithm requires minimum human eflEbrt and little insight 
into the details of chemical mechanism to generate reliable kinetic parameter. 

Based on the above review it is clear that GA has not been applied to ODH 
reaction. Recently, Sulay has used GA to analyze the effect of support monolayer vanadia 
on titania and alumina and obtained promising results. However, the objective fimction 
used in that work is not amendable for multi-response systems. In this study it is observed 
that power law models better explain the data. 

2.3.1 Kinetic studies involving ODH 

Few kinetic studies involving the ODH of propane over metal oxide based catalysts exist. 
A brief review of these studies is given below. 

Andersson [14] studied the partial pressure dependencies of propane oxidation. 
The corresponding rate expression, which is in accordance with Rideal type rate 
expression for various concentrations of vanadium on amorphous AIPO4, is adequate to 
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explain the partial pressure dependencies of both oxygen and propane. A first order 
dependency of propane for the ODH reaction is proposed. 

Chen et al. [16] performed kinetic isotopic studies for the ODH of propane on 
10% V 205 /Zr 02 catalyst. The activation of methylene C-H bond is found to be the 
kinetically relevant step for the ODH of propane, while for secondary combustion allylic 
C-H bond in propene is rate-controlling step. Furthermore, Chen et al [44] proposed that 
the breakings of C-H bond, dissociative adsorption of oxygen are irreversible steps. 

Creaser and Andersson [12] have undertaken a detailed kinetic investigation of 
ODH on V-Mg-0 catalyst. Various rate expressions are derived in order to fit the data 
based on Power Law type models as well as mechanistic models, such as Langmuir- 
Hinshelwood and Mars-Van-Krevelen models. In all the models, they have considered 
carbon oxides as secondary products. However, in a later study Creaser et al. [13] have 
concluded that the reaction data is best explained by considering carbon oxides as 
primary and secondary combustion product. Furthermore, they have proposed 
involvement of lattice oxygen in the ODH reactions. 


2.4 Summary 

From the above review, it is clear that supported vanadium oxide and chromium oxide 
show the promising catalysts for the propane GDH reaction. The literature review of the 
characterization and reaction studies for the supported vanadium and chromium oxide 
catalysts reveal that surface metal oxide species form a unique phase, which is physically 
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and chemically different from the bulk phase. Oxide support and metal oxide loading are 
critical parameters in deciding the nature of the metal oxide species and its reactivity. To 
understand the effect of surface metal oxide phase, CraOs/TiOa and V205/Ti02 are 
considered for propane ODH reaction. Furthermore, to understand the effect of loading, 
several loadings of V205/Ti02 are also studied. The effect of metal oxide phase and metal 
oxide loading is studied with the help of observing the changes, if any, in the kinetic 
parameters. The kinetic parameters are estimated using GA. 
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CHAPTER 3 


EXPERIMENTAL DETAILS 
3.1 Catalyst Preparation 

Catalyst samples used in the present study were prepared by the incipient wetness 
impregnation method. This method is known for its simplicity and viability for making 
low amounts of metal oxide catalysts. Two types of catalysts were prepared; a series of 
supported vanadium oxide (vanadia) catalysts and supported chromium oxide (chromia) 
catalysts. Titania is used as the support. The precursors used for supported vanadia and 
chromia catalysts were ammonium metavanadate (NH 4 VO 3 ) and chromium nitrate 
nanohydrate (Cr(N 03)3 . 9 H 2 O, Aldrich, 99.98% purity), respectively. Prior to preparing 
supported vanadia catalyst the support, titania (Degussa, P-25, 55m^/g), was initially 
pretreated with incipient volumes of known concentration of oxalic acid solution. The 
support was then dried at room temperature in a desiccator followed by drying at 383 K 
for 8 h. The support was finally calcined at 723 K for another 8 h. Prior to making the 
supported chromia catalyst the support was pretreated with incipient volumes of water. 
The heating treatment similar to the above was then followed. 

For the supported vanadia series, the incipient wetness impregnation method 
of the above-pretreated support was carried out with a solution of vanadium oxalate 
precursor. The vanadium oxalate solution was prepared by adding known amounts of 
ammonium meta vanadate with stoichiometric amount of oxalic acid to distilled water. 
Stirring of this mixture was done until a homogeneous solution was formed. The deep 
blue solution formed by the above method was fiirther diluted with double distilled water 
such that the total volume corresponds to the incipient wetness impregnation volume of 
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the support. The pretreated support was then taken in a crucible and the precursor 
solution was added drop wise and mixed. The paste was then heat treated as described for 
the pretreatment of the support. By this method 1, 3 and 4% V2O5 loadings on titania 
support are prepared. Additional details can be found in elsewhere [28]. 

For the chromia-titania catalyst, a predetermined amount of chromium nitrate 
nanohydrate precursor was intimately mixed with incipient wetness volume of water and 
subsequent heat treatment conditions for this catalyst were similar to those described for 
the supported vanadia catalyst. In this way, 1 and 4% Cr203 loadings on titania support 
were prepared. Additional details can be found in elsewhere [36]. For both catalyst x% 
refers to the weight percent of vanadia or chromia, in terms of V2O5 or Cr203, in the 
catalyst. 

3.2 Experimental Set up 

The oxidative dehydrogenation (ODH) of propane was performed in a down flow tubular 
quartz reactor with inlet dia 10 mm and outlet dia 5 mm. The catalyst sample amounts 
used in the reaction ranged from 0.03 to 0.15 g and is given in Table 3.1. Quartz powder 
was used as a diluent and well mixed with supported catalysts in order to prevent 
temperature gradients and to avoid the channeling of gas within the catalyst bed. This 
mixture of catalyst and quartz powder was placed on the well-seated quartz wool at the 
center of the 300 mm long quartz reactor. The catalyst was placed in the isothermal zone 
of the furnace. A schematic view representing the flow diagram of the experimental set 
up is shown in Fig 3.1. 
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The reactant mixture of propane and air at different molar ratios were used with a 
total flow rate of 50 ml/min. Propane to oxygen molar ratios used for both series of 
supported vanadia and chromia samples ranged firom 1:1 to 4:1 and the temperature 
window of operation was from 63 3K to 733K. The flow rate of propane and air was 
controlled by separate mass flow controllers (Bronkhost Hi-Tech, Model F-201D and 
Model F-201C) A Chromel-Alumel thermocouple was inserted from the top of the 
reactor as a sensor in a feedback control loop, central to which was a PID controller (Fuji 
PXZ-4) capable of controlling the reactor temperature within 1 ° C range of accuracy to 
the set point. 

Analysis of product gases was achieved using an on-line gas chromatography 
(NUCON, MODEL-5765), which was equipped with an activated alumina column for 
separation of propane, propene, CO and CO 2 . The GC oven temperature was maintained 
at 327 K and the analysis was performed with the help of a flame ionization detector 
(FED). Carbon oxides, CO and CO 2 , were converted to methane by a methanizer in order 
to be detected by the FED. Main components of the reactor outlet were propane, propene, 
CO and CO 2 . No extra products were observed in case of the vanadia series of samples. 
However, in case of chromia other peaks were observed at high temperatures, which 
could be attributed to oxygenated peaks. Data for each experimental condition was 
collected thrice. Blank runs revealed that there was no propane conversion in the 
experimental region considered suggesting that homogeneous reaction did not take place. 
3.3 Reactivity Calculations 

Analyzed products of the reaction were propene, carbon monoxide, carbon dioxide and 
unconverted propane. Each of these components gave separate areas in the 
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chromatograph. The calculations of conversion, activity, selectivity and yield per gm 
catalyst were based on mole fractions. The area of each component is converted to mole 
fraction based on the response factor given in Table 3.2. 

3.3.1 Conversion 

Conversion for this reaction can be calculated by using following formulae, 

^ ■ n/ Moles of Propane converted 

Conversion, % = x 100 

Moles of Propane in 


Moles of Propane in - Moles of Propane out 
Moles of Propane in 


xlOO 


(3.1) 


3.3.2 Selectivity 


Selectivity for the reactions leading to formation of propene was calculated as follows; 

Propene selectivity, % = .. x loo (3.2) 

Moles of Propane converted 

3.3.3 Yield 


The amount of propene formed was represented by yield, which was calculated as 
follows: 

_ -.JO/ Moles of Propene formed 

Propene yield, % = x 100 (3 .3) 

Moles of Propane in 

3.3.4 Carbon Balance 

Carbon balance represents the accuracy of measurements and provides an indicator as to 
whether carbon is deposited on the catalyst. Carbon balance calculations are performed 
based on the formula given below. 

1 . t - 1 0 / Moles of Carbon out 

Carbon balance, % = ^ — x 1 00 (3 .4) 

Moles of Carbon in 
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3.4 Kinetic Modeling of ODH 

Improvement of kinetic modeling for ODH reaction in the absence of well-accepted 
mechanism is still a challenging task. It initially involves, modeling of this complex 
reaction followed by estimation of kinetic parameters. Subsequently insight for pathways 
of different products can be obtained. 

3.4.1 Parameter Estimation 

Earlier attempts to estimate kinetic parameters include linearisation of the rate equation 
so that ordinary least square methods can be applied. However, it is always not possible 
to linearize rate expression of heterogeneous catalytic reactions and often the estimation 
of rate parameters is achieved by nonlinear regression methods. Consequently, practical 
problems due to non-linearity of the models and limitations of traditional optimization 
methods are usually encountered. The objective function based on nonlinear model for 
the kinetic study of experimental data usually contains more than one optimum and there 
by traditional methods are prone to non-global optima solutions. 

The present study involving the estimation of parameters from propane ODH 
reaction makes use of an evolutionary search algorithms called Genetic Algorithm (GA). 
This is a class of algorithm that is mainly dictated by Darwinian principle of survival of 
the fittest [41]. The successful identification of reasonable parameters gives insight about 
catalyst dependent kinetics of the reaction. Procedure for modeling of ODH reaction, 
outlines of GA technique and different reaction schemes for Power Law type models is 
presented here. 
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3.4. 1.1 Modeling of reaction 

Kinetic analysis of the ODH reaction is performed using the integral method. Reaction 
rate equations with coefficients that have Arrhenius dependence on temperature require 
nonlinear procedure to obtain rate constants. Non-linear regression analysis has the 
following advantages over conventional linear regression approach: 

• To find rA, traditional methods require plots of Xa vs. W/Fao- Drawing the 
tangent at various points gives the rate. This is somewhat subjective in nature and 
usually not reproducible. 

• During linearization, taking reciprocals and cross-multiplication in mechanistic 
models distort the error structure because the rA quantity is in the denominator. 
The error in rAis meaningful when rAis small. 

• Conventional methods are mostly applicable for single reaction. For complex 
reaction systems, however, each of the rates is not known. Therefore in these 
cases the integral approach is useful. 

Other advantages in detail are described elsewhere [46] 

The following assumptions are made for better implementation of the reaction 

• The reaction is operated under steady state isothermal conditions. 

• Total numbers of moles of gases remain constant. 

• Gas phase reactions are negligible. 

• Catalyst activity is not affected due to deactivation. 

These assumptions are justified since propane conversions are maintained less 
than 10%, several runs are taken at each experimental condition and no reaction occurs in 
an empty reactor 
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MATERIAL BALANCE; Material balance calculations based on the above assumptions 
for each component, i, and a particular reaction network can be written as 

Iw ~Kl^i Iw+dM' = 0 Z = 1 p (3.1) 

j 

Where, 

Vo = volumetric flow rate of the feed 
Xi = mole fraction of the i*^ component 

njj = stoichiometric coeflBcient of the i*** component for the j*** reaction 
rj = rate of j**' reaction in terms of volume per unit mass of catalyst per unit time 
and is a function of and x; 

Kj = kinetic parameters for the j*** reaction 
W = weight of the catalyst, 
p = number of components 

Simplification of equation 3.1 using above assumptions gives rise to 


dx^ 

'dW 





i = \ .p 


(3.2) 


The set of ordinary differential equations represented by equation 3.2 can be 
solved to obtain the output mole fraction of each component based on the input value of 
mole fraction, Xio and knowledge of rj. The initial mole fraction is changed during 
integration over the entire mass using Runga-Kutta fourth order technique. The specific 
rate of reaction, rj, is dependent on values of kinetic parameter Kj, and some or all mole 
fractions, Xj. For a particular set of input parameters and the functional dependency, q = 
:^(xi, K^, the integration of the above equation will give the output mole fraction. The 
predicted the output mole fraction obtained this way can then be compared with actual 
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output mole fraction. Minimization of the difference between the actual and predicted 
mole fraction forms the basis for changing the parameter values, Kj, by GA. For multi- 
response systems, the objective function is suitably modified and discussed later in 
section 3. 4.2. 3. 

3.4.2 The concept of Genetic Algorithm (GA) 

Genetic algorithms are different from most classical optimization methods. Principles of 
GA techniques are based on natural genetic inheritance and selection of evolution 
process. These are multidirectional search and optimization procedures that combine 
survival of the fittest among a set of string structures with a structured yet randomized 
information. The string structures are modified based on experience. New sets of artificial 
strings are created in every generation using bits and pieces of the fittest of the old 
artificial strings. These GA’s are known for striking the balance between exploration and 
exploitation of given search space to estimate parameters [41], 

3.4.2. 1 Optimal problem formulation using GA 

Problem formulation procedure is necessary to create a mathematical model for the 
present optimal design problem of parameter estimation. GA requires choosing required 
variables, constraints, objective function and variable bounds. A detailed optimal 
problem procedure for general case is described in [45,47]. 


22 



3.4.3.2 Working principle 

Genetic algorithms possess the advantage of not using an initial estimate of parameters 
for its operation. Some of the features of GA used in the present study are given below. 
CODING: The GA variables are coded (binary or real) from lower to upper bound on the 
variable to form a string. The theoretical foundations of genetic algorithm depend on 
binary string representation of solutions and on the notion of schemata and templates 
allowing exploration of similarities among the chromosomes. 

ASSIGNING FITNESS TO SOLUTION: The objective function evaluated at certain 
values of variables is transformed to non-negative values based on the fitness function. 
The fitness or objective function considered is discussed in section 3. 4.2.3. The GA 
searches for solutions that minimize the fitness function. 

REPRODUCTION OPERATION. The primary objective of the reproduction operator is 
to emphasis good solution and eliminates bad solution while keeping the population size 
constant. Using this operator multiple copies of good solution are made and bad solution 
are eliminated based on a probabilistic approach. 

CROSSOVER OPERATOR: Crossover is the mating process that allows information 
exchange from one generation to the next. In crossover, two randomly chosen strings 
exchange their sub-strings to create two new strings. The number of crossover operations 
is controlled by the crossover probability. Pc, which represents the probability that a 
particular pair of chromosomes is selected for crossover. If the crossover is too high then 
highly fit individuals are discarded faster than selection can produce the improvements. 
On the other hand, if it is too low, selection may stop for lack of exploration. Usually Pc 
value lies between 0.5 and 1.0 [41]. 
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MUTATION OPERATOR: The mutation operator changes a 1 to a 0 and vice a versa in a 
binary string to create a new string, which hopefully have better fitness value. The 
mutation operation is controlled by the mutation probability. Pm. The optimum balance 
between convergence and divergence is maintained by a high Pm and low Pc. The values 
of Pm and Pc used in the present study are 0.006 and 0.9 respectively. 

3. 4.2. 3 Objective function 

To obtain the optimized kinetic parameters, the choice of the objective function is 
critical. Various objective functions like minimization of sum of the squares (SOS) of 
deviations between outlet and inlet concentrations are extensively used in the 
optimization problems. However, in the present case, there are several simultaneous 
inlets and outlet concentration values are present, i.e., the system has multiple responses. 
During situations dealing with multi response systems the minimization of following 
determinant as an objective function is shown to be better [48]. 


det 


(Xji^ expt " Xii^ pred) 2 [(Xii^expt~ Xji^perd) S [(Xii^expt Xil^perd) 2 [(Xil,expt Xil^perd) 
(Xi2,expt ~ Xi2,pred)] (Xi3,e3q)t ~ Xy^pred)] (Xi4,expt “ Xj4,pred)] 

s [(Xii ,expt ” Xii^pred) S (Xi2,expt Xi2,perd) 2 [(Xi2,expt Xi2,perd) ^ [(Xi2,expt Xi2,perd) 

(Xi2,expt “ Xi2,pred)] (^i3,expt 5Ci3^pred)] (Xi4,expt Xi4^pred)] 

[(Xji^expt ■" Xji^pred) 2 [(Xi2,expt “ Xi2,perd) [(Xi3^expt Xi3,perd) 21 [(Xi3^expt Xi3^perd) 

(Xi3,expt “ Xi3,pred)] (Xi3,expt Xi3^pred)] (Xi4,expt Xi4^pred)] 

21 [(Xji^expt "■ Xjljpred) 21 [(Xi2,expt Xi2,perd) 21 [(Xis^expt ““ Xi3^perd) 21 [(XH^expt”" Xi4,perd) 
(Xi4,expt **■ Xi4^pred)] (Xi4^expt Xi4^pred)] (Xi4^expt Xi4^pred)] 
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The problem then can be formulated as minimization of the above determinant with 
respect to Kj, which will give the true estimate of kinetic parameters for a particular 
function dependency, rj = ^(xi, Kj). 

3.4.3.4 Variable bounds 

GA requires certain bounds to be set on each parameter. Optimized parameter values 
should lie within the chosen parameter limits. Ideally GA will access global minimum in 
the bounds of the parameters search space in spite of local minima being present in the 
parameter space. 

3.4.4 Reaction schemes 

The kinetic study of propane ODH includes various empirical and mechanistic models. 
These models proposed can be broadly classified as Power law type models and 
mechanistic models. Most of the models suggest that propene formation is a result of 
primary reaction. However, the formation of carbon oxides as primary or secondary 
reactions is stiU debatable. Here discussion is confined to Power Law based models and 
details about mechanistic models, such as Langmuir-Hinshelwood model, Rideal type 
model, Mars-van-Krevelen model and model involving pseudo-first-order 
approximations are not considered 

The reaction of propane with oxygen can be presented in the form of six reactions 
steps, ri to rs, as shoAvn with the generalized scheme in Fig.3.2. Propene is formed by the 
ODH of propane (ri), which can degrade to form CO and CO 2 by reaction r 2 and rs. 
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respectively. The carbon oxides, CO and CO 2 , can also be formed directly from propane 
by reaction r 4 and rs, respectively. Furthermore, CO 2 can be formed by the oxidation of 
CO by reaction re. Carbon containing compoimds are only ;^ 0 ’wn in Fig 3 .2. 



Fig 3.2 Generalized scheme for Power Law model reaction network 
The stoichiometric equations corresponding to each reaction in the reaction network 
shown in Fig. 3.2 can be written as 


C 3 F 3 +0.50, -+C,H, +H,0 

(3.3) 

+30, -^3C0 + 3H,0 

(3.4) 

+A.50, -+300, +3H,0 

(3.5) 

C,H^+3.50,-^3C0 + AH,0 

(3.6) 

+50,^3 CO, + 4H,0 

(3.7) 

CO +0.50, -+CO, 

(3.8) 


26 





3.4.5 Power law models 


Several power law models were analyzed to determine the empirical reaction rates, which 
can further be used to develop mechanistic models. In the power law models it is 
assumed that the reaction rate is proportional to the partial pressure of reactant raised to a 
particular exponent. Based on the reaction scheme given in Fig. 3.2 the rate expression 
are expressed as 


''2 

^3 


^4 

r. 


^C3HS ^02 


^ 'pni ' 

^C3H6 ^02 


^ ^ jym ' 

^C3H6 ^02 


4 pm ^ ^ pm ^ 

^C3HZ ^02 


,Pn. 


^ ^ pm ^ ^ ^ 

^C3m ^02 


^6 (“)"'(—) 

^ ^ pm '' ^ptw^ 

^CO ^ 02 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


Where, 

ki is the rate constant for reaction i, and 

a; and hi are the partial pressure exponent of the reactant for reaction i 


In these reactions the partial pressures are centered around a particular mean pressure pi “ 


and discussed later. 
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3.5. 1.1 PL-1 model 


The PL-1 model assumes that CO and CO 2 are secondary products and are formed from 
propene. Thus n, r 2 and rs are the only reactions considered since u, rj and re are 
insignificant. This reaction model is shown in Fig 3.3. 



Fig 3.3 PL-1 reaction model 

Based on the reaction scheme given in Fig 3.3 the rate expression are expressed as 


^ ^ ^ jnm ' ^ T>m ' 

^02 


^ -i V pm ' ^ pm ' 


C2E^ ^ 02 


^ ■* /o 


C3H6 ^ 02 


(3.15) 

(3.16) 

(3.17) 
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Rate expression represented by ri, r 2 and ra are the rates of formation of propene, carbon 
monoxide and carbon dioxide respectively. 

3.5.1.2 PL-2 model 

The PL-1 model given above is modified, considering that CO and CO 2 can be primary 
and /or secondary products. Propane reacts to form CO, CO 2 and C 3 H 6 . In this model re is 
insignificant. The reaction scheme is shown in the Fig 3.4. 

CO 



Fig 3.4 PL-2 reaction model 
The rate expression are expressed as 


1 1 v pm ' ^ pm ' 

^02 


r, 

2 2 \ rtm ' ^ niM ' 


Pr 


C2H6 ^ 02 


(3.18) 


(3.19) 
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(3.20) 


■3 ^ ^ pfii ^ ^ pm ' 

^C3H6 ^02 


r - Ir r ^CZHZ \aA r^02 \M 
^ C3HZ ^02 


5 D \ pm ' ^ pm ' 


C3H8 ^ 02 


3.4.5.3 PL-3 model 


(3.21) 

(3.22) 


In this model, it is assumed that CO 2 is formed from CO only. The series reaction 
involved is propane giving rise to propene, which further reacts to form CO. The CO 
formed oxidizes to CO 2 . Thus, ri, ra and re are the only reactions occurring and the rest 
are insignificant. The reaction network is given in Fig 3.5 and corresponding rate 
equations are shown below. 



Fig 3.5 PL-3 reaction model 
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(3.23) 
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(3.24) 


^ ^ pm ^ pm ^ 

^C3H6 ^02 


h=kA^yu^f^ 

o o \ pm / ^ pm 

I nr\ A ^ 


CO ^02 


3.4.5.4 PL-4 model 


(3.25) 


The PL-3 model is further modified to consider a parallel path for CO production from 
propane. The reactions occurring are n, X 2 , u and u and are shown in Fig 3.6. The 
reactions rs and rs are insignificant. The corresponding rate equations for Xi, X 2 , X 4 and re 
are given below. 



Fig 3.6 PL-4 reaction model 
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(3.27) 
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(3.28) 


(3.29) 


n=kA—Y\—T 

O O \ T\m ^ ^ T^m ^ 


CO 


02 


(3.30) 


The number of parameters associated with the various power law models is given in 
Table 3.3. Various parameters for Power Law models and their ranges of search for the 
present discussion of ODH reaction are reported in the Table 3.4. Number of kinetic 
parameters depends on the type of the model selected as well as assumptions made in that 
proposed model. 


3.6 Model selection 

When replicate experiments are not available and in the absence of error analysis, F-test, 
a statistical criterion for simultaneous reactions, is used to select suitable model that best 
explains the reaction data. 
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where. 


yii,= Experimental concentration of h"’ response in i‘'’ experiment 

= Predicted concentration of k'" response in r’ experiment 

p = Number of kinetic parameters in a model 
n = Number of experiments 
v=Number of responses 

cr“ =Elements of inverse of the (v x v) error covariance matrix. (=1) 

Among the set of available models, one with highest Fc will be considered to be the 
best among the set of rival models [48], 
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Fig 3 . 1 REACTOR SETUP FOR PROPANE ODH REACTION 
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Table 3.1: Different weights of catalyst samples used during reaction studies 


Catalyst 

Weight (g) 

1% V205/Ti02 

0.100 

3% V205/Ti02 

0.040 

4% V205/Ti02 

0.030 

l%Cr203/Ti02 

0.150 

4% CrzOa/ Ti02 

0.050 













Table 3.2: Response factor values used in carbon balance. 


Component 

Response Factor 

Propane 

0.6162 

Propene 

0.8744 

Carbon Monoxide 

1.0000 

Carbon Dioxide 

1.0000 
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Table 3.3 Different Power Law models and corresponding number of kinetic parameters 


KINETIC 

MODEL 

CONTRIBURTING 

RATES 

NUMBER OF PARAMETERS 

kfo 

Ei 

ai 

bi 

Total 

PL-1 

ri , T2, rs. 

3 

3 

3 

3 

12 

PL-2 

^ I7 ^2,'i ^ 3 ^ ^45 ^ 5 . 

5 

5 

5 

5 

20 

PL-3 

ri, ra, re. 

3 

3 

3 

3 

12 

PL-4 

ri, ra, T 4 , re. 

4 

4 

4 

4 

16 
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Table 3.4 Lower and upper bounds of the kinetic parameters for Power Law models 


Kinetic 

parameter 

unit 

Bounds 

Lower 

Upper 

Kio. 


0.01 

20 

Ei. 

kjmol ^ 

20 

200 

ai, bi. 

Dimensionless 

0.01 

1.00 
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Chapter 4 

RESULTS AND DISCUSSION 

4.1 Characterization studies 

Characterization studies provide information about the nature of active species. 
Combining spectroscopic studies with reaction data is very handy in order to interpret 
experimental and modeling results. 

The surface areas of V205/Ti02 samples for different loadings are almost constant 
(~ 40 m^/g) for 1 to 4% V2O5 loadings. BET areas for different V205/Ti02 loadings are 
reported in Table 4.1. The XRD studies does not detect separate peak corresponding 
V2O5 phase. However, from Raman spectroscopic studies obtained under ambient 
conditions, the presence of V2O5 crystals are detected for the loadings above ~ 4% 
V205/Ti02 sample [28]. Thus, monolayer coverage is achieved for 4% V205/Ti02. Similar 
kind of Characterization analysis for Cr203/Ti02 samples is performed. BET areas of 
different Cr203/Ti02 loadings are ranged from 40 to 60 mVg (pure titania). These BET 
areas for different Cr203/Ti02 loadings are also reported in the Table 4.1. From Raman 
and other spectroscopic studies, it is observed that monolayer coverage is achieved at ~ 
4% Cr203/Ti02, In situ Raman and DRS characterization studies under dehydrated 
conditions revealed the presence of dispersed surface Cr^^ oxide species for the samples 
below monolayer loadings. Additional Cr^^ and 5 - Cr^^ species were detected by EPR 
[36]. 

Below the monolayer coverage both vanadia and Chromia surface species are 
present in the respective catalysts. The nature of these surface metal oxide species have 
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been studied by various techniques and structures of the surface vanadia and Chromia 
have been proposed [49], 

4.2 Reaction studies 

Reaction studies were carried out over 1, 3 and 4% V 205 /Ti 02 and 1 and 4% Cr 203 /Ti 02 
catalysts. During the reaction studies the total feed flow rate was maintained at 50 cc/min 
and propane to oxygen gas ratio was varied fi-om 1:1 to 2:1 to 3:1 to 4:1. A known 
amount of catalyst was placed in the isothermal zone of the furnace. All the runs were 
taken in the temperature range of 633 to 733 K with 20 K difference between the 
consecutive runs. At each reaction condition data is collected thrice. Based on the inlet 
propane (CsHg) and oxygen (O 2 ) concentrations and outlet areas detected by GC, the 
input and output mole percents of CsHg, CsHg, CO and CO 2 are tabulated in the Tables 
4.2 to 4.6. In these Tables the first column corresponds to reaction temperature. The 
second and third major columns give the input (CsHg and O 2 ) and output (CsHg, CsHe, 
CO 2 and CO) mole percents. The last column gives the C-balance in terms of the ratio of 
carbon atoms out to carbon atoms in, Cout/Cm, and is a measure of the accuracy of 
analysis. The data in each Table is organized such that propane to oxygen ratio of 1:1 is 
the top set of data and is separated by a gap fi'om the following ratio 2:1 data, which is 
followed by the 3:1 and then the 4:1 data. In each data set, corresponding to a particular 
propane to oxygen ratio, the data is presented in increasing reaction temperatures. 

4.2.1 The effect of vanadia loading on activity and propene selectivity 

Previous studies on the same catalysts have suggested that the activity during the propane 
ODH increases with metal oxide loading up to monolayer coverage [28]. To confirm this 
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objective, the activity and propene selectivity for 1, 3 and 4% V205/Ti02, from the data 
given in the Tables 4.2 to 4.4 and are shown in the Figs 4.1 to 4.4. In the Figs. 4.1 to 4.4, 
the molar ratio of C3H8:02 changes from l.T to 2:1 to 3:1 to 4:1. For all these C3H8:02 
ratios, the activity increases with loading for the temperature range used in the present 
study. The catalytic activity also increases with increase in temperature, which is 
expected. Furthermore, for a particular catalyst, the activity appears to increasing with 
C3H8:02 ratio (see, for example 673 K). The propene selectivity, however, appears to 
reach a maximum for the 3% V205/Ti02, except for the C3H8:02 ratio of l.T at 713 K (Fig 
4.1). The propene selectivity decreases with increase in temperature. With a change in 
C3H8:02 ratio, no clear trend in propene selectivity is observed for a particular catalyst. 

4.2.2 The effect of vanadia and chromia species on propane ODH 

The effect of changing the surface metal oxide species on Ti02 was tested by comparing 
the activity and selectivity trends for 1% V205/Ti02 and 1% Cr203/Ti02 catalysts as a 
function of temperature. The data from Tables 4.2 and 4.5 are used to plot Fig 4.5. The 
activity for both catalysts increases with temperature. This increase in activity with 
temperature is expected. The variation in the selectivity with increase in temperature is, 
however, different for the two catalysts. For 1% V205/Ti02 the propene selectivity 
decreases with increase in temperature, whereas, for 1% Cr203/Ti02 the propene 
selectivity increases with the temperature. Thus, surface metal oxide species has a crucial 
role in determining the yield of propene. Consequently, it is worthwhile to study the 
kinetics of the reaction to better understand the factors governing the activity and propene 
selectivity as a function of loading and surface metal oxide species. 
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4.3 Kinetic investigation 

The absence of a known rate determining step and dealing with a set of rival models that 
appear to almost equally fit the experimental data, kinetic investigation, therefore, is a 
twofold problem of parameter estimation and model discrimination. Ideally, the estimated 
parameters, which are statistically significant, should be positive and comply with 
physicochemical constraints. In the present study the best model selection is based on the 
statistical F-test criteria for simultaneous reactions [48]. 

4.3.1 Parameter estimation 

As discussed previously, the estimation of kinetic parameters for oxidative 
dehydrogenation of propane is done using genetic algorithms. Parameter estimation is 
based on the power law type models, whose detailed description is made in experimental 
section 3.4. The five catalysts considered were 1, 3 and 4% VsOs/TiOa and 1 and 4% 
Cr203/Ti02. 

In order to study the effect of loading for V205/Ti02 catalyst, parameter 
estimation for various power law models is done using the data from Tables 4.2 to 4.4 for 
the 1, 3 and 4% V205/Ti02 catalysts, respectively. Kinetic parameter estimation is based 
on the minimization of multiresponse objective function for simultaneous reactions. For 
parameter estimation, the inverse temperature is centered about the mean reaction 
temperature and its value is observed to be 683.16 K. This centering is applied to the 
concentration terms in the form of Pm- The Pm values for propane and oxygen are 0.328 
and 0.141 atm, respectively. These Pm values for the propene and carbon monoxide are 
reported in the Tables containing the kinetic parameters. Kinetic parameter estimation is 


42 



done for all models of the present study and the parameters for each catalyst are 
presented in the Tables A, 1 to A.4 for the four power law models. 

4.3.2 Model selection 

The set of simultaneous rate equations, dictated by the type of proposed model, 
determines the kinetic behavior of the ODH reaction scheme. When replicate experiments 
are not available and at the absence of pure error analysis, the criteria for selection of the 
best model that best explains the reaction data among the set of rival models for 
simultaneous reactions can be done using F-test. 

Selection of a suitable model or models is dictated by the critical value. The F“ 
values are obtained from the published data and depend on the significance level, number 
of parameters, number of variables and number of experimental runs. Based on the data 
and the model, the Fo value can be calculated. For a particular catalyst, a model is 
statistically significant provided, the Fo value is greater than the F*^ value. Furthermore, 
among the Fo values obtained for different models, the one with the highest Fo value is the 
most suited model to describe the experimental data. The Fo and F“ values are given in the 
Tables 4.7 and 4.8 for V205/Ti02 and Cr203/Ti02 catalysts, respectively. From the Table 
4.7 it is observed that all the models are statistically inadequate to describe the 
experimental data for the 1% V205/Ti02 catalyst since the Fo value is less than the F°o.io 
(F'" value at 10% significant level). Though, none of the models is statistically adequate, 
the highest value of Fo is obtained for the PL-1 model. In this model carbon oxides are 
formed as secondary products. For the 3 and 4% V205/Ti02 catalysts it is observed that 
Fo values for all the models are greater than F'^o.os values. In fact, the Fo values are greater 
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than F‘'o.oi values for both the catalysts, except for the PL -2 model of 4 % V205/Ti02 
catalyst. Similar to the 1 /o V 2 O 5 /Ti 02 catalyst the experimental data obtained for the 3 % 
V205/Ti02 is best represented by the PL -1 model. For 4 % V205/Ti02 catalyst, however, 
the PL -4 model best represents the experimental data. In the PL -4 model, CO is a primary 
and secondary product and CO2 is produced only from CO. 

The Fo and values for the Cr203/Ti02 catalysts are presented in the Table 4 . 8 . 
From the Table, it is observed that for the 1 % Cr203/Ti02 catalyst, the Fo values are 
greater than F^o.os values, except for the PL -2 model. Similar to the sub monolayer 
V205/Ti02 catalyst, the Fo value for PL -1 model is highest and, thus, best explains the 
experimental data. For the 4 % Cr203/Ti02 catalyst the Fo values are better than even the 
F'^0.01 value. Thus, all models are statistically significant. Among the different models, the 
highest Fo value is obtained for the PL -3 model, which suggests that this model best 
represents the experimental data for the 4 % Cr203/Ti02 catalyst. 

4.3.3 Modeling analysis based on proposed best models 

Once the most suited models have been determined based on the highest Fo values the 
kinetic analysis can be done in more detail. The analysis is done with respect to the 
loading and effect of surface metal oxide species for the best models. 

The effect of vanadia or chromia loading is first analyzed with respect to ODH 
reaction, ri and the respective parameters kio, Ei, ajand bi. These values are given in the 
Table 4.9 for the four catalysts. The values of kio are observed to increase with loading 
for V205/Ti02 and Cr203/Ti02 catalysts. Specifically, for the best models kio values for 
the 3 and 4 % V205/Ti02 catalyst and the 1 and 4 % Cr203/TiO2 catalysts are 13 . 25 , 14 . 48 , 
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2.34 and 7.46 in the respective units. Values of kio are roughly scale per amount of active 
metal oxide site. Furthermore, comparing the two different metal oxide species it is 
observed that the kio values for V205/Ti02 catalysts are greater than the Cr203/Ti02 
catalysts. Thus, the higher kio values suggest that V205/Ti02 catalysts are suited to 
produce propene. The ranges of activation energies are, however, in different region. For 
the V205/Ti02 catalysts the activation energies are 61 and 58 kJ/mol, whereas, for the 
Cr203/Ti02 catalysts the activation energies are 87 and 131 kJ/mol. From the activation 
energy values of V205/Ti02 catalysts, it can be inferred that same reaction path for 
propane ODH exists both for sub monolayer and monolayer catalysts as there is 
negligible difference in the activation energy values, however, for Cr203/Ti02 catalysts, it 
appears that propane ODH path for sub monolayer is different from monolayer catalyst as 
there is considerable difference in activation energy values (44 kJ/mol). Furthermore, 
propane and oxygen partial pressure dependencies, ai and bi, are different for the four 
catalysts. The propane partial pressure dependency exists for the 3 and 4 % V205/Ti02 
and 1 % Cr203/Ti02, whereas, no propane dependency exists for the 4 % Cr203/Ti02. The 
4 % V205/Ti02 and 4 % Cr203/Ti02 are strongly oxygen dependent, whereas, 1 % 
Cr203/Ti02 and 3 % V205/Ti02 are weakly oxygen partial pressure dependent. 

The information about the degradation of propene can be, at the first 
approximation, related to the values k2o, the pre-exponential factor for CO formation. For 
the best models values of k2o, E2, a2 and b2 are reported in the Table 4 . 10 . For the best 
models, the values of k2o are 9 . 54 , 7 . 02 , 1.76 and 15.14 for 3 and 4 % V205/Ti02 and 1 
and 4 % Cr203/Ti02 catalysts. The pre-exponential factors for the propene degradation are 
1/3 of these values since three moles of CO are produced from the each C3H8. Thus, the 
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difiference between the kio and k 2 o /3 is greater for the V205/Ti02 catalysts compared to 
the Cr203/Ti02 catalysts. Consequent^, the net formation of propene on V205/Ti02 
catalyst is larger compared to the Cv^O-i/TiOi catalysts. The activation energies for the 
propane degradation reaction, r2, range from 50 to 71 kJ/mol for the four catalysts. 
Comparison of Ei and E2 vdues provides information regarding the temperature 
dependency of the related rates of propene formation and propene degradation. The 
difference between the Ei and E2 are larger for the Cr203/Ti02 catalysts than the 
V205/Ti02 catalysts. This suggests that with an increase in temperature the rate of 
formation of the propene would increase for the Cr203/Ti02 catalysts, whereas, the rate 
of formation of propene would increase more gradually or perhaps decrease with an 
increase in temperature. Indeed, the propene selectivity decreases with temperature for 
the V205/Ti02 catalysts and increases with temperature for the Cr203/Ti02 catalysts as 
seen in Fig 4 . 5 . 

Thus, based on the kinetic parameters it is possible to explain the effect of the 
loading and surface metal oxide species. The values of the kio and k2o govern the rate of 
propene formation and the activation energies Ei and E2 determine the effect of the 
propene selectivity with temperature. The difference in the values of Ei between the 
V205/Ti02 and Cr203/Ti02 catalysts suggests that reaction pathways for the important 
ODH reaction are different. 
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Table 4.1: BET area values for V 205 /Ti 02 and CxiOj,/ Ti02 catalysts 


Catalyst 

Surface Area (m^/g) 

Reference 

Ti02 

60 

Ref[28] 

1% V205/Ti02 

39 

Ref [28] 

3% V205/Ti02 

41 

Ref [28] 

4% V205/Ti02 

40 

Ref [36] 

1% CrsOs/ Ti02 

41 

Ref [36] 

4% CrzOs/ Ti02 

n.d 



n.d not determined 
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Table 4.2: Input and output mole percent of components for 1% V 205 /Ti 02 catalyst. 
Total flow rate = 50 ml/min. Weight of the catalyst =100 mg. 


Temperature 

(K) 

Input mole (%) 

Output mole (%) 

C- 

balance 

Cout/Cin 

C 3 H 8 

O 2 

C 3 H 8 

C 3 H^ 

CO 2 

CO 

673 

17.355 

17.355 

17.037 

0.112 

0.164 

0.087 

0.99 

713 

17.355 

17.355 

16.700 

0.307 

0.181 

0.177 

0.98 

733 

17.355 

17.355 

16.650 

0.448 

0.258 

0.293 

0.99 


653 

29.577 

14.789 

29.495 

0.215 

0.103 

0.063 

1.00 

673 

29.577 

14.789 

27.637 

0.256 

0.132 

0.152 

MEM 

713 

29.577 

14.789 

26.299 

0.611 

0.924 

0.404 

WEME 


653 

38.650 

12.883 

38.451 

0.159 

0.089 

0.097 

1.00 

693 

38.650 

12.883 

38.073 

0.330 

0.422 

0.243 

0.99 

713 

38.650 

12.883 

35.550 

0.395 

0.494 

0.363 

0.99 

733 

38.650 

12.883 

34.000 

0.950 

0.597 

0.519 

0.96 


653 

45.652 

11.413 

45.405 

0.217 

0.122 

0.116 

1.00 

673 

45.652 

11.413 

43.892 

0.298 

0.420 

0.201 

0.97 

693 

45.652 

11.413 

42.894 

0.479 

0.715 

0.369 

0.96 















Table 4.3: Input and output mole percent of components for 3% V 205 /Ti 02 catalyst. 
Total flow rate = 50 ml/min. Weight of the catalyst = 40 mg. 


Temperature 

(K) 

Input mole (%) 

Output mole (%) 

C- 

balance 

Cout/Cin 

CsHg 

O 2 

C 3 H 8 

C 3 H 6 

CO 2 

CO 

653 

17.355 

17.355 

16.620 

0.297 

0.132 

0.135 

0.98 

673 

17.355 

17.355 

16.574 

0.383 

0.189 

0.222 

0.98 

693 

17.355 

17.355 

16.399 

0.476 

0.335 

0.361 

0.99 

713 

17.355 

17.355 

15.900 

0.567 

0.385 

0.545 

0.97 


673 

29.577 

14.789 

28.768 

0.501 

0.165 

0.195 

0.99 


29.577 

14.789 

28.240 

0.652 

0.263 

0.332 

0.98 

713 

29.577 

14.789 

27.501 

0.836 

0.499 

0.522 

0.97 


673 

38.650 

12.883 

37.380 

0.543 

0.183 

0.223 

0.98 

693 

38.650 

12.883 

37.353 

0.729 

0.306 

0.381 

0.99 

713 

38.650 

12.883 

36.620 

0.892 

0.556 

0.622 

0.98 


633 

45.652 

11.413 

45.051 

0.483 

0.135 

0.176 

1.00 

673 

45.652 

11.413 

44.291 

0.933 

0.488 

0.610 

1.00 

693 

45.652 

11.413 

43.290 

1.098 

0.771 

0.955 

0.98 

713 

45.652 

11.413 

43.242 

1.366 

1.195 

1.166 

0.99 

733 

45.652 

11.413 

41.877 

1.628 

1.172 

1.520 

0.97 
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Table 4.4: Input and output mole percent of components for 4% V 205 /Ti 02 catalyst. 
Total flow rate = 50 ml/min. Weight of the catalyst = 30 mg. 


Temperature 

(K) 

Input mole (%) 

Output mole (%) 

C- 

balance 

Cout/Cin 

CsHs 

O 2 

C 3 H 8 

C 3 H 6 

CO 2 

CO 

653 

Mmai 

17.355 

17.014 

0.218 

0.071 

0.170 

0.99 

673 

tsm 

17.355 

16.624 

0.423 

0.407 

0.532 

1.00 

713 

17.355 

17.355 

16.143 

0.772 

1.31 

1.137 

1.00 


653 

29.577 

14.789 

28.464 

0.443 

0.120 

0.188 

0.98 

673 

29.577 

14.789 

27.278 

0.655 

0.428 

0.407 

0.95 

713 

29.577 

14.789 

27.072 

0.914 

1.117 

0.723 

0.97 


633 

38.650 

12.883 

38.665 

0.354 

0.169 

0.149 

1.00 

653 

38.650 

12.883 

37.916 

0.649 

0.362 

0.300 

1.00 


38.650 

12.883 

36.708 

0.769 

0.506 

0.479 

0.98 

713 

38.650 

12.883 

35.970 

1.056 

1.140 

1.211 

0.98 


633 

45.652 

11.413 

44.916 

0.472 

0.125 

0.187 

0.99 

653 

45.652 

11.413 

43.655 

0.513 

0.537 

0.390 

0.97 

673 

45.650 

11.413 

42.602 

0.820 

0.816 

0.658 

0.96 

693 

45.652 

11.413 

41.341 

1.134 

1.110 

1.030 

0.95 
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Table 4 . 5 : Input and output mole percent of components for 1 % Cr203/Ti02 catalyst. 
Total flow rate = 50 ml/min. Weight of the catalyst =150 mg. 


Temperature 

(K) 

Input mole (%) 

Output mole (%) 

C- 

balance 

Cout/Cin 

CsHg 

O2 

CsHg 

C3H6 

CO2 

CO 

633 

17.355 

17.355 

16.768 

0.170 

0.125 

0.140 

■DESH 

653 

17.355 

17.355 

16.103 

0.250 

1.732 

0.232 

■iK:M 

673 

17.355 

17.355 

15.543 

0.380 

2.462 

0.383 

Hsai 


673 

29.577 

14.789 

27.994 

0.171 

0.741 

0.130 

0.96 

693 

29.577 

14.789 

27.715 

0.217 

0.984 

0.210 

0.96 

713 

29.577 

14.789 

27.313 

0.299 

1.220 

0.378 

0.95 









653 

38.650 

12.883 

38.441 

0.047 

0.682 

0.321 

1.00 

673 

38.650 

12.883 

37.731 

0.245 

0.813 

0.356 


693 

38.650 

12.883 

37.068 

0.731 



1.00 

713 

38.650 

12.883 

34.741 

1.100 

2.290 

1.039 

0.96 


673 

45.652 

11.413 

45.065 

0.361 

1.256 

0.337 

1.00 

653 

45.652 

11.413 

44.605 

0.496 

1.368 

0.527 

1.00 

673 

45.652 

11.413 

44.319 

0.818 

1.602 

0.835 

0.99 













Table 4.6: Input and output mole percent of components for 4% Cr 203 /Ti 02 catalyst. 
Total flow rate = 50 ml/min. Weight of the catalyst = 50 mg. 


Temperature 

(K) 




c- 

balance 

Cout/Cin 

CsHs 






673 

17.355 

17.355 

15.460 

0.712 

1.66 

0.630 

0.98 


653 

29.577 

14.789 

28.700 

0.348 

0.98 

0.267 

0.99 

693 

29.577 

14.789 

28.260 

0.760 

1.310 

0.540 

WESM 

723 

29.577 

14.789 


1.511 

1.990 


1.03 


623 

38.65 

12.883 

38.490 

0.043 



1.00 

653 

38.650 

12.883 

37.770 

0.335 

0.922 

0.246 

0.99 


38.-650 

12.883 

37.018 

0.420 

0.938 

0.304 

0.98 

693 


12.883 

37.300 

0.530 



0.99 

723 

38.650 

12.883 

34.430 

1.200 

1.400 

0.700 

0.96 


653 

45.652 

11.413 

44.530 

0.311 

1.296 

0.257 

0.99 

673 

45.652 

11.413 

44.060 

0.378 

1.370 

0.340 

0.98 

693 

45.652 

11.413 

43.739 

0.629 

1.500 

0.462 

0.98 





















Table 4.7: Modeling (Fo) and critical (F“) values of F-test for titania supported 
vanadium oxide catalysts in propane ODH reaction esc] . 


Model 

1% V205/Ti02 

3% V205/Ti02 

4% V205/Ti02 

B 





Fo 

F^O.05 

F^.oi 

PL-1 




IQH 

2.496 

2.57 

1.930 

2.520 

PL-2 

1.11 




2.223 

2.14 

1.770 

2.247 

PL-3 

1.42 

1.675 

n 

3.25 

1.917 

2.496 

2.78 

1.930 

2.520 

PL-4 

1.13 

1.607 

3.06 

1.815 

2.315 

2.94 

1.828 

2.339 


F^*o.o 5 = Critical value ofF at 5% significant level. 
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Table 4.8: Modeling (Fo) and critical (F*^) values of F-test for titania supported 
chromium oxide catalysts in propane ODH reaction 


Model 

1% CrzOa/TiOa 

4% Cr203/Ti02 

Fo 

F^O.05 

F*^ 0.01 

Fo 

F%.05 

F'^o.oi 

PL-1 

3.05 

1.944 

2.547 

4.168 

1.960 

2.579 

PL-2 

1.22 

1.790 

2.275 

3.01 

1.807 

2.307 

PL-3 

2.61 

1.944 

2.547 

4.65 

1.960 

2.579 

PL-4 

1.93 

1.842 

2.366 

3.65 

1.859 

2.399 











Table 4. 10: BCinetic parameters from the best suited models for the degradation of 
propene to carbon monoxide formation 


VzOs/TiOa 

KINETIC 

PARAMETER 

Cr203/Ti02 

1% 

3% 

4% 

1 % 

4% 

0.75 

9.54 

7.02 

k 20 ml STP min^ (g cat)'^ atm ' 

1.76 

15.14 

61 

50 

60 

E 2 kJ mol^ 

67 

71 

0.16 

0.66 

0.01 

CI 2 

0.3 

0.28 

0.01 

0.01 

0.99 

b2 

0.01 

0.91 

PL-1 

PL-1 

PL-4 

Best suited model 

PL-1 

PL-3 
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Figure 4.1; Variation of activity and propene selectivity with vanadia loading for propane 
ODH reaction at different temperatures. Total flow rate = 50 ml/min. Molar 
ratio of CsHg: O 2 = 1:1 
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Figure 4.2; Variation of activity and propene selectivity with vanadia loading for propane 
ODH reaction at different temperatures. Total flow rate = 50 ml/min. Molar 
ratio of CsHg : O 2 = 2:1 
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Figure 4.3: Variation of activity and propene selectivity with vanadia loading for propane 
ODH reaction at different temperatures. Total flow rate = 50 ml/min. Molar 
ratio of C3H8:02=3:1 
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Figure 4.4: Variation of activity and propene selectivity with vanadia loading for propane 
ODH reaction at different temperatures. Total flow rate = 50 ml/min. Molar 
ratio of CsHg : O 2 = 4:1 


fin 
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Figure 4.5: Variation of activity and propene selectivity 'with temperature for 1% V205/Ti02 
and 1% Cr20.i/Ti02 in the propane ODH reaction. Total flow rate = 50 ml/min. 
Molar ratio of C3H8;02 = 3:1. 

Legend: O- Propene selectmty for 1% V205/Ti02 catalyst, •-Activity for 1% 

V205/Ti02 catalyst, O -Propene selectivity for 1% Cr203/Ti02 catalyst, 
Actmty for 1% Cr203/Ti02 catalyst. 


Propene selectivity (%) 



Chapter 5 


CONCLUSIONS AND RECOMMENDATIONS 

Estimation of kinetic parameters using GA in order to explain the effect of surface metal 
oxide species (vanadium oxide and chromium oxide) on titania support gave the insight 
about the propane ODH reaction kinetics. 

5.1 Conclusions 

Based on the reaction and modeling studies in the present study, the following 
conclusions can be drawn. 

1 . From the reaction studies, it is observed that propane conversion increase with 
loading for V205/Ti02 and Cr203/Ti02 in the region considered. This can be 
attributed to increase in number of two-dimensional active species, which are 
favorable to ODH of propane to propene. The effect of temperature, however, is 
different. For Cr203/Ti02 catalysts propene selectivity increases with temperaure 
and for V205/Ti02 catalysts propene selectivity decreases with temperature. 

2 . Using genetic algorithm the kinetic parameters for different power law models are 
successfully obtained. 

3 . For both vanadia and chromia catalysts in the ODH of propane, it appears sub 
monolayer loadings follow the PL -1 model and monolayer loadings follow PL -3 
or PL -4 models. 

4 . The best suited models suggest that monolayer V205/Ti02 and Cr203/Ti02 
catalysts favour the oxidation of carbon monoxide to carbon dioxide 

5 . Considering the activation energy for ODH of propane to propene, for V205/Ti02 
catalyst, the same path for ODH reaction exists for sub monolayer and monolayer 



catalysts. However, for CraOs/TiOa catalyst, propane ODH reaction path changes 
from sub monolayer ( 1 % Cr203/Ti02) to monolayer ( 4 % Cr203/Ti02) catalyst. 

6 . Based on the kinetic parameter values the effect of loading on the propane activity 
and effect f the temperature on the propene selecthity can be explained. 


5.2 Recommendations 

The following recommendations are proposed for ftirther study; 

1 . Further work on the effect of loading and nature of the surface metal oxide 
species should be extended to mechanistic models, which provide more 
information about the kinetics. 

2 . Detailed information can be extracted using advanced in situ techniques to find 
the active sites responsible for both activity and selectivity. Correlate the activity 
and selectivity of these catalysts with in situ characterization techniques may give 
more realistic information. 

3 . Efficient modeling can be achieved in the ODH reaction by considering the effect 
of water CO, CO2 and C3H6 in the inlet. 
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appendix 1 


Reference factor values are 


obtained using following formula 


For hydrocarbons, the 


relative FID response ftctor is stated as. 


^ aw 


)/ u xH „ ) 
^ nw n 


X 0.7487 


Where; 

F, = relative response factor for a hydrocarbon 


or hydrocarbon group type 


of the same 


carbon number. 

Ctrw= atomic weight of carbon, 12.011- 

c,= number of carbonmolecules in thehydrooarbon. 


atomic eight of hydrogen, 1 . 008 , 

number of hydrogen molecules in the hydrocarbon, 
0.7487 - corrects the response of methane to umty. 
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rable A. 1 : Kinetic parameters corresponding to PL-1 model for different catalysts 


KINETIC 

V205/Ti02 

Cr203/Ti02 

PARAMETERS 

1% 

3% 

4% 

1% 

4% 

kio 

1.59 

13.25 

18.87 

2.34 

13.17 

k20 

0.75 

9.54 

10.52 

1.76 

6.92 

ho 

1.18 

10.84 

10.33 

5.27 

15.05 

El 

79 

62 

66 

87 

77 

E2 

61 

50 

64 

67 

42 

Es 

24 

41 

86 

41 

25 

ai 

0.58 

0.58 

0.51 

0.72 

0.01 

bi 

0.52 

0.01 

0.29 

0.13 

0.99 

a2 \ 

0.16 

0.66 

0.45 

0.30 

0.41 

b2 

0.01 

0.01 

0.01 

0.01 

0.99 

as 

0.35 

0.88 

0.41 

0.14 

0.19 

bs 

0.99 

0.26 

0.46 

0.97 

0.03 

Fitness Value 

6.69E-021 

2.86E-021 

4.02E-20 

4.34E-019 

8.02E-022 

Pc3H6 

3.67E-3 

7.50E-3 

4.40E-3 

4.30E-3 

6.00E-3 

Generation 

1000 

1000 

1000 

1000 

1000 

Mass (g) 

0.1 

0.04 

0.03 

0.15 

0.05 
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able A.z: Kinetic parameters corresponding to PL-2 model for different catalysts 


KINETIC 

PARAMETERS 

VzOs/TiOz 

Cr203/Ti02 

1% 

3% 

4% 

1 % 

4% 

kio 

1.23 

13.01 

16.50 

2.26 

6.97 

k 20 

0.10 

7.33 

8.61 

0.50 

1.30 

kso 

0.06 

6.36 

4.75 

2.04 

4.74 

k40 

0.62 

2.42 

0.70 

1.63 

2.95 

kso 

1.00 

4.41 

3.90 

4.74 

10.20 

El 

85 

59 

56 

101 

115 

E 2 

142 

57 

70 

20 

105 

Es 

SI 

29 

20 

86 

29 

E 4 

62 

69 

193 

100 

82 

Es 

52 

95 

165 

77 

45 

aj 

0.59 

0.59 

0.67 

0.67 

0.29 

bi 

0.03 

0.01 

0.56 

0.56 

0.99 

ao 

0.01 

0.53 

0.07 

0.82 

0.01 

bz 

0.99 

0.01 

0.01 

0.29 

1 0.99 

as 

0.13 

0.99 

0.57 

0.09 

0.90 

bs 

0.99 

0.45 

0.01 

099 

0.71 

a4 

0.99 

0.99 

0.99 

0.94 

0.01 

bs 

0.01 

0.01 

0.01 

0.99 

0.54 

as 

0.98 

0.99 

0.04 

0.01 

0.99 

bs 

0.60 

0.73 

0.96 

0.99 

0.81 

Fitness Value 

1.24E-021 

2.08E-021 

2.46E-020 

3.67E-019' 

8.01E-022 

Pc3H6 

3.67E-3 

7.50E-3 

4.40E-3 

4.30E-3 

6.00E-3 

Generation 

1000 

1000 

1000 

1000 

1000 

Mass (g) 

0.1 

0.04 

0.03 

0.15 

0.05 


69 




Table A.3: Kinetic 


inetic parameters corresponding to PL-3 model for different catalysts 
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Table A.4: Kinetic parameters corresponding to PL-4 model for different catalysts 


KINETIC 

PARAMETERS 

1% 

VaOs/TiOa 

3% 

VsOs/TiOa 

4% 

V 2 O 5 / TiOs 

1% 

Cr203/Ti02 

4% 

Cr203/Ti02 

kio 

1.66 

11.96 

14.49 

3.06 

7.13 

^20 

3.36 

10.77 

7.02 

8.62 

4.15 

k40 

0.37 

6.46 

6.49 


16.35 

k60 

2.19 

11.61 

9.17 

4.79 

15.46 

El 

66 

56 

58 

69 

106 

E 2 

20 

49 

60 

34 

25 

K 4 

20 

89 

176 

20 

67 

E/5 

20 

20 

97 

56 

20 

ai 

0.19 

0.67 

0.82 

0.66 

0.01 

bi 

0.01 

0.06 

0.99 

0.01 

0.99 

a 2 

0.63 

0.61 

0.01 

0.51 

0.60 

b2 

0.02 

0.32 

0.99 

0.59 

0.74 

04 

0.39 

0.99 

0.17 

0.28 

0.01 

b4 , 

0.12 

0.01 

0.01 


0.97 

a<5 

■ 0.37 

0.80 

0.35 

0.46 

0.44 

be . 

0.49 

0.12 

0.50 

0.24 

0.15 

Fitness Value 

2.32E-021 

2.58E-021 

3.43E-020 

4.63E-019 

1.55E-021 

Pc3H6 

0.0036 

0.0075 

0.0044 

0.0043 

0.0060 

Pco 

0.0024 

0.0053 

0.0054 

0.0042 

0.0044 

Generation 

1000 

1000 

1000 

1000 

1000 

Mass (g) 

0.1 

0.04 

0.03 

0.15 

0.05 
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